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ABSTRACT 


This  paper  introduces  a  new  tool  for  image  texture  analysis 
called  a  polarogram.  A  polarogram  is  a  polar  plot  of  an  orientation 
sensitive  texture  statistic.  Polarograms  give  rise  to  a  class  of 
texture  descriptors  which  are  sensitive  to  both  texture  coarseness  and 
directionality,  but  yet  which  are  invariant  to  rotations  of  the  image 
textures.  An  experiment  is  described  in  which  polarograms  are  applied 
to  the  classification  of  texture  samples. 


Accession  For 

ATI'S  GRAM 
DDC  TAB 
Unannounced 
Justification 


Distribution/ _ _ 

Availability  Codes 

Avail  and/or 
Dist.  special 


$ 


DTIC 

ELECTE 
JUN  10  1980 


■■  isms.zn  (apsc) 

I'cvU.ved  and  is 

l^J-18  (7b). 


ko,  ;.:s  v  •  •  .  '  '  .  ; 

ii. t  tf. 

c  't'.'  „ 

Di-  \.rii .  . .  ‘.id. 

A.  i), 

l.-ehnlool  Infora&Vlaa  Officer 


1  •  introduction 

This  paper  introduces  a  new  tool  for  image  texture 
analysis  called  a  polarogr am.  A  polarogram  is  a  polar 
representation  of  an  orientation  specific  texture  statistic. 
It  can  be  used  to  obtain  information  about  both  texture 
coarseness  and  texture  directionality. 

Texture  coarseness,  or  the  size  of  the  texture  elements, 
is  ordinarily  related  to  the  autocorrellation  of  the  texture. 
For  example,  the  contrast  feature  for  grey  level  cooccurrence 
matrices  (Haralick [ 1] )  or  grey  level  difference  histograms 
(Weszka  et  al[2])  is,  essentially,  a  measurement  of  a  single 
value  of  the  image  autocorrelation.  It  is  ordinarily  compu  :ed 
for  pairs  of  image  points  which  are  close  to  one  another  in 
the  image.  High  values  of  contrast  usually  indicate  fine 
textures  (where  there  are  relatively  many  edge  points)  and  low 
values  of  contrast  indicate  coarse  textures.  Of  course,  other 
factors  such  as  edge  sharpness  also  effect  the  value  of  the 
contrast  statistic. 

Texture  directionality  has  not  been  handled  as  adequately 
as  coarseness.  It  has  been  suggested  that  directionally 
specific  statistics  be  computed.  For  example,  suppose  that 
P(o,r)  is  a  polar  representation  of  the  power  spectrum  of  a 
texture.  Then  Weszka  et  a 1 ( 2 J  suggested  "wedges"  defined  as: 


w(n1,n2) 


as  direction  sensitive  texture  statistics.  A  difficulty  with 
such  statistics  is  that  they  are  not  rotation  invariant  - 
i.e.,  rotating  the  original  texture  changes  the  statistics. 
Since  orientation  cannot,  in  general,  be  controlled,  and  might 
even  vary  across  a  single  "homogeneous"  field,  it  is  important 
that  even  statistics  which  measure  directionality  be  rotation 
invariant.  Polarograms  can  be  used  to  generate  rotation 
invariant  statistics  which  are  sensitive  to  texture 
directionality. 

In  Section  2  we  define  the  polarogram,  and  present  some 
examples.  Section  3  contains  an  experimental  study  using  the 


texture  database  employed  in  [3j.  Finally,  Section  4  contains 
conclusions  and  a  discussion. 


. . .  . 
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2 . POLAROGRAMS 


Polarograms  are  a  new  tool  for  describing  image  textures. 
A  polarogram  is  a  polar  plot  of  a  texture  statistic  as  a 


function  of 

orientation.  For 

example , 

let 

D(a) 

be  a 

displacement 

vector  of  fixed 

magnitude , 

d , 

and 

var iable 

orientation,  a,  and  CD^aj  be  the  cooccurrence  matrix  for 
displacement  D(a).  Let  f  be  some  statistic,  such  as  contrast, 
defined  for  a  cooccurrence  matrix.  Then  we  can  define  the 
polarogram,  P^,  by 

Pf(a)  =  f(CD(a)) 

As  a  simple  example,  consider  the  vertical  bar  texture  in 
Figure  la.  The  bars  are  five  pixels  wide.  Let  D(a)  have 
magnitude  1,  and  let  f  be  the  contrast  statistic.  Then  Figure 
lb  displays  the  polarogram  Pf.  Note  that  Pf  (  tt  /2)  =  0  since 
all  pairs  of  adjacent  vertical  points  have  identical  grey 
levels.  P f  ( 0 )  ,  on  the  other  hand,  has  value  .4,  since  2  out 
of  every  5  one's  are  adjacent  to  a  zero.  In  general,  if  we 
regard  the  texture  in  Figure  la  as  a  continuous  image,  then 
Pf(a)  =  2cos(a)/5.  In  practice,  Pf(a)  is  only  computed  for  a 
discrete  set  of  values  for  a.  An  interpolation  function  is 
used  to  compute  P^(a)  for  intermediate  values  of  a.  The 
simplest  interpolation  scheme  connects  consecutive  values  of 
Pj(a)  with  straight  lines. 

Texture  statistics  are  derived  from  polarograms  by 
computing  size  and  shape  features  of  the  polarogram.  The 
shape  features  will  ordinarily  not  only  depend  on  the  shape  of 
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the  boundary  of  the  polarogram,  but  also  on  the  position  of 
the  origin  of  the  polarogram  within  that  shape.  For  example, 
the  polarograms  in  Figures  2a  and  2b  have  the  same  boundary 
shapes  (circles) ,  but  the  polarogram  in  Figure  2a  is  centered 
at  the  circle  center,  while  the  polarogram  in  Figure  2b  is 
centered  away  from  the  circle  center.  This  difference 
indicates  a  crucial  difference  in  the  two  underlying  textures. 
The  texture  corresponding  to  the  polarogram  in  Figure  2a  is 
isotropic  with  respect  to  the  statistic  of  the  polarogram, 
while  the  texture  corresponding  to  Figure  2b  is  not.  The 
skewness  feature  defined  below  measures  the  extent  of  such 
isotropy. 

*\n  important  class  of  statistics  which  can  be  computed 
from  a  polarogram  are  its  moments  from  the  origin  (note  that 
these  are  not  the  same  as  the  central  moments  of  the 
polarogram).  The  p'th  moment  of  P^  is  : 


P 


f 


(a)  da 


i 


where 
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In  the  special  case  where  the  boundary  of  P^  is  a 
polygon,  it  can  be  shown  that: 

U1  =  al 

2 

u2  “  a2  “  U1 

3 

U3  —  a 3  3u^ aj  ^  2u^ 


where 


and 


yi/2’apli) 


i=l 


a1(i) 


a2(i) 


a3(i) 


a^b^  siny^ 


c.  y . 

1  x 


log 


(ci+ai  -  b.cosYi)  (ci4bi  -  a.cosYi) 

.  .  2 
a^b^sxn  y^ 


.  /  a.b.  .  \  2 

=  - —  (  1  —  sln^iJ  (cot  a.  +  oot  B-) 

V  C.  /  1  1 

i  i a,b'  \ 3  r 

=  (  -i— -  siny.  ]  cosec  a.  cot  a. 

2y.V  c.  x)  [  xi 

+  log  |  |(cosec  +  cot  ou )  (cosec  B^  +  cot  B^) 


+  cosec  8^  cot  Bj 


See  Figure  3  for  the  definitions  of  n b  ^  y^,  a^,  b^,  c^. 

The  first  moment,  u^,  represents  the  deviation  from  the 
mean  of  the  perimeter  distance  to  the  origin,  U2  the  variance, 
and  U3  the  skewness  of  the  perimeter  distribution.  Two  simple 
measures  of  size,  which  are  used  in  the  experiments  in  Section 
3,  are  area  of  the  polarogram  and  perimeter. 
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Notice  that  the  statistics  computed  from  the  polarogram 
are  invariant  to  the  orientation  of  the  polarogram  -  i.e.,  if 
P'£<a)  =  Pf(a+da),  then  P  and  P'  will  have  the  same 
statistics.  Since  rotating  the  original  texture  results  in 
rotating  the  polarograms,  these  statistics  are  invariant  to 
rotations  of  the  image  textures,  yet  can  still  be  used  to 
measure  texture  "directionality." 

As  an  example,  Figure  4  contains  two  texture  samples,  one 
of  grating  and  the  other  of  metal  scrap.  Figure  5  contains 
polarograms  for  the  textures  in  Figure  4.  These  are  based  on 
the  contrast  statistic  for  cooccurrence  matrices  using  D(a) 
with  fixed  magnitude  s  e  symmetry  in  the  polarograms  is 
due  to  the  symmetry  in  the  cooccurence  matrices  -  CD^a)  = 
CD(a+pi).  Cooccur rence  matrices  were  computed  for  directions 
which  are  multiples  of  45  degrees.  The  larger  size  of  the 
grating  polarogram  reflects  the  higher  contrast  of  the  grating 
(note  that  the  first  order  grey  level  statistics  of  the  two 
textures  are  identical) .  The  relatively  low  value  of  the 
grating  polarogram  in  the  horizontal  direction  reflects  the 
elongation  of  the  grating  texture  elements,  an  aspect  of 
texture  directionality.  The  scrap  metal,  on  the  other  hand, 
is  more  isotropic,  so  its  polarogram  is  more  circular. 


i 


Figure  4b.  A  64x64  metal  t 


Figure  5b.  Polarogram  for  ! 
and  the  contrast  descriptor 


a)  grating  (G) 


b)  concrete  (E) 


c)  pebbles  (P) 


d)  tree  bark  (T) 


e)  metal  (M) 


P  f  i  rn  f) . 


T nvf-}iro  q.omnl  r<; . 
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3.  EXPERIMENTS 

An  experiment  was  performed  to  illustrate  the  use  of 
polarogram  statistics  for  texture  classification.  A  database 
of  five  texture  classes  was  chosen,  and  ten  64x64  samples  were 
included  in  each  class.  The  textures  in  the  database  include 
grating  (G) ,  tree  bark  (T) ,  metal  scrap  (M) ,  pebbles  (P) ,  and 
concrete  (E) .  These  samples  form  part  of  a  larger  database 
described  in  experiments  reported  in  [3].  Figure  6  contains 
one  64x64  sample  from  each  of  the  five  classes.  All  of  the 
textures  have  been  subjected  to  a  grey  scale  normalization 
procedures  which  "flattens"  their  histograms  and  guarantees 
that  all  of  the  samples  have  identical  first-order  grey  level 
statistics . 

Polarograms  were  computed  for  the  contrast  descriptor  of 
the  grey  level  cooccurrence  matrix.  Cooccurrence  matrices 
were  computed  for  each  of  the  eight  principle  directions  on 
the  digital  grid,  and  for  distances  between  pairs  of  points  of 
1,  3  and  5.  Thus,  three  polarograms  were  computed  for  each 
texture  sample.  For  each  of  the  polarograms,  the  five 
statistics  discussed  in  Section  2  were  computed.  Three 
classification  experiments  were  performed;  one  based  on  the 
distance  1  statistics,  one  based  on  the  distance  3  statistics 
and  one  based  on  the  distance  5  statistics. 

The  textures  were  classified  using  a  multivariate  linear 
discriminant  function  [4],  computed  based  on  the  assumptions 
that  the  probability  density  functions  for  the  statistics  of 


all  classes  were  multivariate  normal,  with  a  common  covariance 
matrix,  but  different  mean  vectors,  and  that  the  prior 
probability  of  each  class  is  known.  Subroutine  ODNORM  of 
IMSLIB  was  use  to  compute  the  discriminant  function  and  to 
perform  the  classification. 

The  results  of  the  experiments  are  summarized  on  Tables 
1-3.  Both  the  distance  5  and  the  distance  3  polarograms  gave 
classif ication  rates  of  88%  on  this  database  of  50  textures. 


G 

E 

P 

T 

M 


G 

1.0 

0 

0 

.  1 
0 


E 

0 

1.0 

0 

0 

0 


p 

0 

0 

.9 

.1 

.2 


T 

0 

0 

0 

.7 

0 


M 

0 

0 

.1 

.1 

.8 


Percent  error  =  12% 


Table  2.  Confusion  matrix  for,  D=3 
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4.  SUMMARY 

We  have  presented  a  new  computational  tool  for  image 
texture  analysis  called  a  polarogram,  and  illustrated  its  use 
by  applying  it  to  a  texture  classification  problem. 
Polarograms  were  designed  to  yield  texture  descriptors  which 
are  sensitive  to  both  texture  coarseness  and  directionality. 
An  important  aspect  of  directionally  sensitive  polarogram 
statistics  is  that  they  are  invariant  to  rotations  of  the 
image  texture.  We  are  currently  applying  the  polarogram  to 
the  unsupervised  segmentation  of  scenes  containing  many 
differently  textured  areas. 
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